Potassium ground state scattering parameters and 
Born-Oppenheimer potentials from molecular spectroscopy 
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We present precision measurements with MHz uncertainty of the energy gap between asymptotic 
and well bound levels in the electronic ground state X 1 E^" of the 39 K2 molecule. The molecules are 
prepared in a highly collimated particle beam and are interrogated in a A-type excitation scheme 
of optical transitions to long range levels close to the asymptote of the ground state, using the 
electronically excited state A Ej as intermediate one. The transition frequencies are measured 
either by comparison with 1 2 lines or by absolute measurements using a fs-frequency comb. The 
determined level energies were used together with Feshbach resonances from cold collisions of 39 K 
and 40 K reported from other authors to fit new ground state potentials. Precise scattering lengths 
are determined and tests of the validity of the Born-Oppenheimer approximation for the description 
of cold collisions at this level of precision are performed. 

PACS numbers: 32.10.Fn, 31.30.Gs, 33.20.Kf 
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I. INTRODUCTION 

Progress in experiments with cold, ultracold, quantum 
degenerate atomic or molecular ensembles starting from 
laser cooled atoms is steady and fast. An appropriate 
model of binary collisions is essential for many current 
and future experiments. The scattering length charac- 
terizes the phase shift of the atom pair wave function 
when one atom slowly propagates through the potential 
induced by the other atom. This parameter reduces the 
interatomic potential to an effective phase shift in the 
limit of zero kinetic energy. It can describe cold colli- 
sions, e.g. the dynamics of thermalization or the stability 
of Bose-Einstein condensates, and is a helpful parameter 
in the understanding of phase diagrams of diluted ultra- 
cold ensembles. 

However, a parametrization of the scattering process 
only by the scattering length is often not sufficient, when 
tunable resonances are involved that e.g. alter the scat- 
tering length [l| or allow the formation of quantum de- 
generate ensembles of dimers Q and demand thorough 
investigations of phase transitions [3, |j] . In many cases, 
models for ultracold collisions were obtained from reso- 
nance spectroscopy even without deep knowledge of the 
interatomic potentials at small internuclear separations 
SI- 

Measuring the binding energy of few loosely bound 
molecular levels may yield enough information to derive 
the scattering length [7J] or Feshbach resonances, but one 
relies on long-range interaction parameters calculated 
from atomic properties. It is also possible to measure 
the energy position of a larger number of weakly bound 
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molecular levels with molecular spectroscopy, i.e., rela- 
tive to a deeply bound molecular level @, Q , or by pho- 
toassociation spectroscopy relative to the atomic asymp- 
tote [nj, LUl QjJ • With such an approach, the determi- 
nation of the asymptotic energy with respect to deeply 
bound levels and the long-range interaction coefficients of 
the potential are reliable if the number of observed levels 
is high and their binding energy is small. This evaluation 
can be combined with calculated long-range parameters 
from dispersion theory since even with a multitude of 
observed states the values of long range parameters are 
significantly correlated and thus difficult to determine ex- 
perimentally in a reliable way. 

If all available data are combined, ground state po- 
tential energy curves can be constructed for the whole 
interatomic separation axis from the repulsive range via 
the bound one to the asymptote, and they contain the 
full information of the molecular level structure and can 
be used to calculate scattering parameters and resonance 
positions. With the power of the full-potential method 
being developed in a preceding experiment on sodium 
[131 ] . the work presented here aims in searching for lim- 
its of the potential method. The Born-Oppenheimer ap- 
proximation (BOA) underlies the potential description 
for molecules. In this approximation, different isotopes 
are described by the same potentials but by the use of 
different reduced masses, hyperfine and Zeeman parame- 
ters in the Hamiltonian. Therefore, it is possible to scale 
the energy levels, Feshbach resonance positions, and scat- 
tering lengths from one isotope combination to another. 
However, this assumes the BOA and any deviation after 
the corrections for hyperfine structure and mass differ- 
ence has to be considered a violation of it. 

Potassium has not only the two stable bosonic isotopes 
39 K and 41 K but also the long-living, naturally occur- 
ring 40 K. The large number of accessible isotope combi- 



nations makes potassium a good candidate for investiga- 
tions of BOA violations. 40 K is of great interest due to 
its fermionic character [g, LLJ] . Feshbach resonances have 
been measured to a high precision for this isotope in sev- 
eral studies [H, [H, O M, GJ, and recently also for 
39 K [21(. Moreover, photoassociation spectroscopy for 
39 K was reported by Wang et al. [Il(, and a rich molecu- 
lar data set obtained mostly on 39 K2 is available for both 
the singlet [22| and the triplet ground state [23| . A small 
number of vibrational levels close to the asymptote could 
not be observed, which makes the global description by 
the existing data incomplete and prohibits conclusions 
about violations of the BOA. 

We report here the observation of levels within the un- 
observed energy interval in a multi-photon process. In 
our study, we use long-range levels of the first electron- 
ically excited state A 1 S+ that we recently investigated 
[241 ] and for which we found indications of non-zero cor- 
rections to the BOA for levels very close to the asymptote 
of this state [2||. We would like to point out that BOA 
corrections discussed therein and in this work are not due 
to hyperfine coupling between different levels, since these 
perturbations are accounted for in the applied Hamilto- 
nian. Non-adiabatic corrections of the BOA are observed 
for chemical reactions of light reaction partners are re- 
ported for spin-orbit coupling by Li Che et al. [2y]. Es- 
timates of adiabatic corrections of the BOA, in a similar 
way like the normal mass effect of atomic energy levels, 
indicate for the molecular potential energy curves that 
extended models for cold collisions might be needed at 
a level of precision of 0.1G or below for Feshbach res- 
onances and molecular levels spanning an energy range 
of several thousand wave numbers with a precision of 
0.001 cm -1 or below as presented here. 

The laser-spectroscopic experiments of this study are 
described in Sec. [TTJ followed by a detailed analysis of 
the frequency measurement in Scc. lIIII In Sec. lIVI we in- 
troduce our model with parameterized potential energy 
curves, which allows us to calculate rovibrational levels 
including hyperfine structure and Zccman effect, Fesh- 
bach resonance positions and scattering lengths. The 
already existing and newly obtained experimental data 
are used to derive quantitative results of the potential 
curves, as discussed in Sec. [V] Moreover, we investigate 
to which precision the BOA is valid by comparing the 
resulting potentials from fits of different combinations 
of isotopic data, thus the mass dependence of the inter- 
atomic interaction potentials. 



II. EXPERIMENTAL SETUP 

The setup, which we use to create a beam of K2 
molecules, was described in some detail in a previous 
publication [271 ] and was well characterized for precision 
measurements [28l . 12911 and Ramsey-Borde interferome- 
try with molecules [30]. In short, potassium is heated 
in a furnace to about 700 K. Through a 200 /im nozzle, 
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FIG. 1: (Color online) Simplified excitation scheme of K2. Ini- 
tially, the molecules are in the vibrational ground state of the 
X J E+ state. The laser Li is applied first and drives Franck- 
Condon pumping to populate v — 41, J in the ground state. 
Asymptotic levels are observed as dark resonances induced by 
laser L2 (on resonance) and laser L3 (tuned). 



potassium vapour expands into the vacuum and forms an 
atomic and molecular beam with a mean particle velocity 
of about 1000 m/s. 

Mainly, the vibrational ground state t>x = of the 
electronic state X 1 T,^ is populated. Due to very un- 
favourable Franck-Condon factors, no direct coupling 
scheme in a A-configuration to asymptotic levels of the 
electronic ground state can be driven starting from this 
level. Therefore, we employ optical pumping in a first in- 
teraction zone of the molecular beam. Laser Li transfers 
a significant amount of population to higher vibrational 
levels [1J, [24], [25| of the ground state. This pumping step 
is achieved by an extended cavity diode laser (ECDL) at 
792 nm in Littrow configuration. The light is transfered 
with a fibre to the vacuum system, where it excites a 
rotational line P(Jx) or R( Jx) in the 23-0 band of the 
A-X system of 39 K 2 (see Fig.[T]). Spontaneous decay pop- 
ulates with about 15% probability the vibrational level 
v = 41, J = Jx-, Jx — 2 or Jx, Jx + 2, respectively. 

Molecules in excited vibrational levels like v = 41 have 
a very long lifetime and will reach the second interac- 
tion zone about 30 cm downstream the molecular beam, 
where they interact with two co-propagating laser beams 
L2 and L3, which are also delivered by fibres. Laser 
fields and the molecular beam interact under a 90° an- 
gle. The laser beams have diameters of a about 2 mm 
and are superimposed on a dichroic mirror. These two 
fields will establish the A-type coupling scheme connect- 
ing the ground state vibrational levels v and v" via the 
vibrational level v' in the excited state A ^^ The en- 
ergetic positions of the vibrational levels in the A state 
were investigated previously and are tabulated in [241 ]. 

Laser induced fluorescence can be observed from the 
second interaction zone. The laser L2 is held on resonance 
for one of the molecular transitions (1/, J' = J ± 1) <— 



(v, J) at about 710 nm and induces fluorescence, mainly 
at 775 nm. L 2 is a titanium:sapphirc laser pumped by a 
10 W pump laser at 532 nm. It is frequency stabilized 
by an offset-lock to an iodine-stabilized HeNe laser. For 
this purpose, L 2 is locked to a tunable transfer cavity, 
which is itself stabilized to the HeNe laser. An acousto- 
optical modulator (AOM) in the beam path of the HeNe 
laser allows a shift of one transverse mode of the transfer 
cavity to any arbitrary frequency of L 2 . 

Laser L3 is a second Tksapphire laser with a 5 W pump 
laser. It is operated around 775 nm and tuned. If L3 is 
resonant with a transition (v 1 , J') — ► (v" , J" = J' ± 1), 
a dark resonance or electromagnetically induced trans- 
parency is induced and the fluorescence from the excited 
state is reduced. Unfortunately, the fluorescence spec- 
trally overlaps with stray light of laser L 3 . Thus, dis- 
crimination by optical filters is not possible and, there- 
fore, direct observation of the dark resonances was not 
successful. We could overcome this problem of detecting 
the resonances in presence of laser stray light by modulat- 
ing laser Li with a chopper and filtering the fluorescence 
signal with a lock-in amplifier driven at the chopper fre- 
quency. 

The signal-to-noise ratio was still small, so the laser L3 
had to be tuned over the dark resonance with high reso- 
lution and very good reproducibility to be able to aver- 
age several independent measurements for improving the 
signal-to-noise ratio. Additionally, its frequency must be 
precisely known (as for laser L 2 ) to calculate the energy 
difference between different asymptotic levels v", J" and 
the levels v, J. We met these requirements by means of 
a two step offset stabilization of L3 with respect to an 
iodine stabilized diode laser (see below) of well known 
frequency. The frequency of laser L 2 was measured in a 
later experimental run by a fs-frequency comb referenced 
to a GPS-disciplincd quartz oscillator. A similar scheme 
as for L3 was not possible for L 2 because no I 2 lines are 
known with the required accuracy in this spectral inter- 
val. 

Figure shows the offset-stabilization scheme of L 3 . 
Iodine lines of the (0-14) and (1-14) bands of the B-X 
system with an uncertainty of about 1 MHz were used 
as references. Though the I 2 spectrum is very dense, 
frequency intervals of tens of GHz had to be bridged to 
reach the resonances of L3. For this purpose we stabi- 
lized an ECDL L rc f by saturated absorption spectroscopy 
to the required transitions. We used 70 cm long iodine 
cells heated to 740 K with a side arm temperature of 
295 K. In principle, L3 could have been stabilized and 
tuned by stabilizing the beat note of lasers L3 and L re f 
to an appropriate tunable high frequency source. How- 
ever, the latter was not available. Instead, we introduced 
an ECDL transfer laser Lt, which is closer in frequency 
than 1 GHz to L3. Lt is frequency stabilized like L 2 
via another transfer cavity relative to the I 2 -stabilizcd 
HeNe laser. The beat note between L ro f and Lt was ob- 
served with a fast photodiode, counted by a rf-counter, 
and recorded for later corrections. The beat note be- 
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FIG. 2: (color online) Laser system used for frequency mea- 
surement and tuning of L3. OD: optical diode; BS: beam 
splitter; PBS: polarizing beam splitter; PZT: piezo-electric 
actuator; A/4: A/4 retarder; AOM: acousto-optical modula- 
tor. 



tween L3 and Lt was mixed with a tunable synthesizer 
and the mixing product was stabilized by means of a 
digital divider and a phase-locked-loop. In this setup we 
transfer the long-term stability and accuracy from the I 2 
lines to the laser L3 and achieve reliable tuning of the 
laser. This laser tuning scheme has already proved to be 
successful for high precision frequency measurement with 
potassium atoms |28| . 



III. MEASUREMENTS AND UNCERTAINTIES 

A. Data Recording 

The aim of the present spectroscopic measurements 
was to determine energy differences between asymptotic 
levels of 39 K 2 and differences of asymptotic levels to levels 
in the well bound region of the potential. For this pur- 
pose, the investigations were divided into two steps: the 
spectroscopy of dark resonances with frequency determi- 
nation of L3 while L 2 was held on a molecular resonance. 
Later we measured with a fs-frequency comb the abso- 
lute frequency of those transitions that turned out to be 
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FIG. 3: (Color online) Spectra of dark resonances observed in two A-configurations. Left-hand side: coupling between (v = 
41, J = 7)-(v' = 143, J' =8)-(u" = 80, J" = 7). Right-hand side: case (v = 41, J = 7)-(w' = 159, J' = 6)-(«" = 83, J" = 5), 
two dips are visible, which belong to the small hyperfine splitting of asymptotic levels in the state X. The hyperfine structure 
stems from coupling to the a 3 Ej state and is so small in potassium that it is only resolved for v" — 82 and higher. The 
frequency scale is given with respect to iodine lines 



successful for the spectroscopy of the dark resonances. 
Systematic frequency shifts between the two sets of mea- 
surements have to be avoided. The main source of un- 
certainty is here the first-order Doppler effect introduced 
by an angle different from 90° between the laser beams 
and the molecular beam (see Sec. IIII B|) . 

We started searching for ground state resonances in 
the region of v" = 80 where the measurements from clas- 
sical spectroscopy [22J stopped and left an energy gap of 
about 1.8 cm" 1 to the asymptote 4s -f 4s. Due to selec- 
tion rules, for each level (v' , J') two dark resonances with 
J" = J' ± 1 are expected for each vibrational level v" . In 
most cases, it is practically impossible to observe more 
than one vibrational level v" from a given level v' be- 
cause the Franck-Condon overlap drops off dramatically 
with the rapidly growing classical outer turning point of 
the vibration wave function when v" is increased by one 
unit. 

In Fig. [21 observed dark resonances for two levels 
(v" , J") are shown. A single resonance appears (left-hand 
side) for the coupling between (v = 41, J = 7) — (v' = 
143, J' = 8) - (v" = 80, J" = 7) when L 3 is tuned. On 
the right-hand side, two resonances are visible that are 
due to hyperfine splitting of the level (v" = 83, J" = 5). 
The structure in the asymptotic levels of the state X 1 E+ 
appears by hyperfine coupling to the a 3 E+ state and is 
small for K 2 compared to Na 2 [M] or Cs2 [31|- We could 
only resolve the splitting for vibrational levels v" = 82 
and higher. 

The hyperfine component at lower frequency in the 



right-hand part of Fig. [3] corresponds to a total nuclear 
spin 1 — 3, the high frequency one to I = 1. Since 
J" is odd, I = 0, 2 do not appear for symmetry reasons 
[3(J | . The intensity ratio of the two observed dips in the 
spectrum is very close to the expected ratio of 7/3. 

We have have observed in total 14 asymptotic levels of 
which the least bound one is only 0.182 cm" 1 below the 
lowest hyperfine dissociation threshold and has an outer 
turning point of about 21 A. The assignment of quan- 
tum numbers to all levels (v",J") is based on selection 
rules for the rotation, and counting of vibrational levels 
provided by Amiot [22j, |32|, which was later confirmed 
by our analysis (Sec. [Vj) . Frequencies, uncertainties and 
quantum numbers of all resonances measured by tuning 
laser L3 are given in Tab. Q] 

After the determination of the transition frequencies 
observed by L3 we have measured the frequency of L2, 
when exciting the upper level (v' , J'), with the help of a 
commercial fs-fibre laser frequency comb, which was ref- 
erenced to a GPS-disciplined quartz oscillator. To avoid 
systematic frequency shifts from first-order Doppler ef- 
fect, laser L2 was aligned with the molecular beam by 
the same procedure (Sec. IIII B|) as before and was sta- 
bilized with the offset-lock to the iodine-stabilized HeNe 
laser to stay resonant with the molecular transition under 
investigation. 

Each frequency setting of laser L2 was typically mea- 
sured for 80 s with a gate time of the counters of 1 s. The 
observed scatter of the measured frequencies was below 
50 kHz. The mode number of the frequency comb was 



TABLE I: Quantum numbers and frequencies v of transitions 
(V — v") observed by tuning L3. If the hyperfme splitting 
could be resolved, the quantum number of the total nuclear 
spin I is given. The number in brackets behind the frequency 
v is the uncertainty in units of the last digit. In the last 
column, the iodine reference lines are listed [33| . 
,, : „-i 



Line 



Iodine line 



143-80 


R(5) 




12 813.819 53(7 


P(72) 0-14 ai 


143-80 


R(7) 




12 813.929 19(7 


P(72) 0-14 ai 


143-80 


R(9) 




12 814.062 36(7 


P(72) 0-14 ai 


143-80 


P(ll) 




12 813.76165(7 


P(72) 0-14 ai 


147-81 


R(l) 




12 843.431 73(7 


R(42) 0-14 ai 5 


147-81 


P(3) 




12 843.371 59(7 


R(42) 0-14 ai 5 


147-81 


R(5) 




12 843.587 27(7 


R(42) 0-14 ai S 


147-81 


P(7) 




12 843.433 71(7 


R(42) 0-14 ais 


147-81 


R(9) 




12 843.845 56(7 


R(42) 0-14 ais 


147-81 


P(ll) 




12 843.605 99(7 


R(42) 0-14 ai 5 


159-82 


R(9) I = 


3 


12 910.763 55(8 


R(102) 1-14 ai 


159-82 


R(9) I = 


1 


12 910.763 84(9 


R(102) 1-14 ai 


159-82 


P(ll)/ = 


= 3 


12 910.587 20(8 


R(102) 1-14 ai 


159-82 


P(ll)/ = 


= 1 


12 910.587 52(9 


R(102) 1-14 ai 


159-82 


P(7)/ = 


3 


12 910.486 21(8 


R(102) 1-14 ai 


159-82 


P(7)J = 


1 


12 910.486 41(9 


R(102) 1-14 ai 


159-83 


R(5) / = 


3 


12 910.100 56(8 


R(102) 1-14 ai 


159-83 


R(5) I = 


1 


12 910.10133(8 


R(102) 1-14 ai 



TABLE II: Quantum numbers and frequencies v of the tran- 
sitions pumped by L2. The number in brackets behind the 
frequency v is the uncertainty in units of the last digit. 



v — V 


Line 


v in cm i 


143-41 


P(7) 


14 057.684 43(7) 


143-41 


R(7) 


14 057.984 44(7) 


143-41 


P(ll) 


14 055.057 67(7) 


147-41 


P(3) 


14 090.628 40(7) 


147-41 


P(7) 


14 089.036 59(7) 


147-41 


R(9) 


14 088.166 78(7) 


159-41 


R(9) 


14 155.964 91(7) 


159-41 


P(7) 


14 157.018 56(7) 



unambiguously determined by measuring the wavelength 
of L2 with a commercial wavelength meter corresponding 
to 30 MHz accuracy. The wavemeter was calibrated by 
measuring a laser stabilized to the Mg intercombination 
line 3J| or to the iodine stabilized HeNe laser used earlier 
for L3. 

The frequencies, uncertainties and quantum numbers 
of all transitions pumped by L2 are given in Tab.[ll] From 
the frequencies given in Tab. U and Tab. [TO various energy 
intervals within the manifold of rovibrational levels of the 
ground state X 1 S+ can be calculated that were used for 
the fits yielding potentials with high precision (Sec.lV|). 



the adjustment of the laser beams L2 and L3 under a 
90° angle to the molecular beam. Any deviation from 
perpendicular adjustment will introduce systematic shifts 
due to first-order Doppler effect. 

We always followed the same alignment procedure in 
each series of measurements. Laser L2 was tuned to the 
potassium Dl line. The atomic line was recorded at low 
laser intensity to avoid saturation broadening and in two 
distinct configurations: with a single pass laser beam and 
with the laser beam retro-reflected. In the second config- 
uration, the beam was reflected back by a mirror about 
2 m from the fibre output. The parallelism of the two 
counter-propagating beams was checked with an aper- 
ture in front of the fibre coupler. 

Then the observed line width in both configurations 
was compared. A deviation from 90° in the crossing an- 
gle of laser and particle beam will show up in a broad- 
ening of the line profile in the two-beam configuration. 
Systematic variation of the angle allowed to determine a 
minimum of the line broadening with an uncertainty of 
about 0.5 mrad, which corresponds to an uncertainty in 
the frequency of L2 of 0.7 MHz due to the Doppler shift. 

For the alignment of L3 with respect to L2, an aperture 
in front of the mirror retro-reflecting L2 was used. Again, 
an uncertainty of 0.5 mrad in the angle was achieved, 
which leads to an additional possible Doppler shift of 
0.6 MHz after alignment of L3 relative to L2. 

The second important contribution to the uncertainty 
is the reliable determination of the resonance frequency 
for laser L2 because the frequency measurements of L2 
and L 3 were separated in time and the dark resonance 
is observed at the two-photon resonance. We found that 
we can determine the line centre to within 1.8 MHz by 
tuning the AOM frequency of the offset-lock of L2 to the 
HeNe laser. This uncertainty is much larger than the 
observed drift of the applied offset-lock during this or 
previous experiments [28j |. 

Further contributions stem from the determination of 
the line centre of the dark resonance (typically 0.3 MHz, 
depending on the signal-to- noise ratio), the frequency un- 
certainty of the iodine lines used for calibration (0.1 MHz, 
for the 1-14 band 1.3 MHz), and the possible uncertain- 
ties in the lock of the spectroscopy laser to the trans- 
fer laser, synchronisation of the counter traces with the 
spectra etc. (in total 0.3 MHz). The uncertainty of the 
frequency measurement with the fs-comb was estimated 
to 0.2 MHz. Other factors like Zeeman effect or second- 
order Doppler effect are negligible. The respective uncer- 
tainty contributions for the measured lines were added 
quadratically. A compilation of all contributions is listed 
in Tab. |Hll 



B. Uncertainty Considerations 



IV. THEORETICAL DESCRIPTION 



Several aspects have to be considered to determine the 
uncertainty of the measured transition frequencies. The 
largest contribution is expected from the uncertainty of 



To cover the full range of experimental data from 
Fourier transform spectroscopy over molecular beam 
spectroscopy to Feshbach resonances by a theoretical 



TABLE III: Summary of the different contributions to the 
uncertainty budget. 



source 


uncertainty 


Doppler effect L2 


0.7 MHz 


Doppler effect L2 rel. to L3 


0.6 MHz 


centre of excitation L2 


1.8 MHz 


centre of dark resonance L3 


0.3 MHz 


offset lock L3, counter 


0.3 MHz 


iodine lines 


0.1 MHz / 1.3 MHz 


fs-comb 


0.2 MHz 



model, only a coupled channel analysis is adequate, which 
includes both molecular states X *E+ and a 3 S+, the hy- 
pcrfinc coupling and the Zeeman interaction, to compute 
the rotational, vibrational and hyperfme structure of the 
bound states and the scattering states. Such a theoret- 
ical approach has been described in several papers, e.g. 
[35|] , so we invite the reader to follow our construction 
of the theoretical model that is able to reproduce all ex- 
isting measurements within their reported experimental 
uncertainty. 

The full Hamiltonian, see for example Refs. [35l. [3rl [371] . 
for a pair of atoms A and B with electron and nuclear 
spin s and i, respectively, can be written in the form: 



H = T n + U x {R)Px + U a (R)(l - P x ) 

+a A (R)s A iA + a B (R)s B iB 

+ {g s AS z A - giAizA)^BohrB z + {g s BS z B ~ 5i-B*z-B)A*Boii 

+~\(R)(3S% -S 2 ) 

The first line shows the kinetic energy T n and the po- 
tential energy Ux and U a for the motion of the atoms, 
where Px and 1 — Px are projection operators on the 
uncoupled states X and a, respectively. The second line 
shows the hyperfme interaction, determined mainly by 
the Fermi contact term, and the functional dependence 
with R for the hyperfme parameters indicates that distor- 
tions of the considered atom from the second one should 
be expected. But the coupling of nucleus A with the 
electron spin of atom B is neglected, and also the inter- 
action with the nuclear quadrupole moment coming into 
play for more deeply bound levels. The third line gives 
the Zeeman energy from the electron spin and the nu- 
clear spin by an external homogeneous magnetic field B 
in z direction. The last line contains the spin-spin in- 
teraction represented by the total molecular spin S and 
its direction to the molecule fixed axis Z. The parameter 
A is a function of R, mainly 1/i? 3 as dipole-dipole in- 
teraction, but it contains also contributions from second 
order spin-orbit interactions. The final analysis showed 
that such later contribution is insignificant within the 
achieved experimental accuracy in the case of potassium. 

We take the magnetic hyperfme parameters and the 
electronic and nuclear g-factors for the atomic ground 
state of the potassium isotopes from the report in [381 ] . 



because the hyperfme structure was only seen for asymp- 
totic vibrational levels and Feshbach resonances, thus any 
dependence of the hyperfme interaction on the internu- 
clear separation by chemical bonds are probably very 
small and not identified in our analysis. The atomic 
masses are taken from the recent tables by G. Audi et 

al. [ai]. 

The calculations incorporate Born-Oppenheimer po- 
tentials according to Hund's coupling case b, because the 
total electronic spin is taken as good quantum number 
for Ux and U a . 

Each representation of the potentials for the two states 
with the common atomic asymptote 4s + 4s is split 
into three regions on the internuclear separation axis R: 
the repulsive wall (R < R ma ), the asymptotic region 
(R > R ou t), and the intermediate region in between. The 
analytic form of the potentials in the intermediate range 
is described by a finite power expansion with a nonlinear 
variable function £ of internuclear separation R: 



Z(R) 



R — Rm 

R + bR r , 



C/i R (i?) = ^a^(P) J , 



(2) 



(3) 



j=0 



where the ai are fitting parameters and b and R m are 
chosen such that only few parameters a,i are needed for 
describing the steep slope at the short internuclear sep- 
aration side and the much smaller slope at the large R 
sjijte by the analytic form of Eq. $5§ , R m is normally close 
to the value of the equilibrium separation. The potential 
is extrapolated for R < R[ nn with: 



U SR (R) = A + B/R 



N a 



(4) 



by adjusting the A and B parameters to get a continuous 
transition at i?i nn ; the final fitting uses N s equal to 12 and 
6 for X : S+ and a 3 E+ states, respectively, as adequate 
exponents. 

For large internuclear distances (R > i? ut) we adopted 
the standard long range form of molecular potentials: 



U LR {R) = U oc - C 6 /R 6 - C 8 /R 8 - C w /R w ± £ eX ch, (5) 
where the exchange contribution is given by 



Scxch = A ox iT exp(-0R) 



(6) 



and Uoo set to zero for fixing the energy reference of the 
total potential scheme. 

The data on hand include three different isotopomers, 
namely 39 K 2 , 40 K 2 , and 39 K 41 K. Thus it is possible to 
check the validity of the Born-Oppenheimer approxima- 
tion or the so called mass scaling of cold collision proper- 
ties between different isotopomers. For this purpose the 



mechanical Schrodinger equation for the rovibrational 
motion should be reformulated to allow for adiabatic and 
non-adiabatic corrections to the Born-Oppenheimer ap- 
proximation. The most often used form (see e.g. [33J) of 
the radial equation for a rotational state J is: 



H 



cir 



+ 



h 2 ■ [l + a{R)\ • J(J + 1) 



2^R 2 



(J) 



For a homonuclear molecule the corrections a, (3 and 
Cad depend on the ratio m e /n of the electron mass m e 
and the reduced mass \i of the molecule. Because only 
two of these corrections can be determined simultane- 
ously [4(|, we decided to neglect in our approach. In 
this way a(R) becomes an effective non-adiabatic cor- 
rection, not only the rotational correction. Additionally, 
the Schrodinger equations gets the proper form for di- 
rect application of the conventional Numerov procedure 
to calculate eigenvalues and eigenfunctions. The adia- 
batic correction C/ a d(i?) can only be determined, if data 
of different isotopic species of the molecule are included. 
Both a and U & d are functions of the internuclear separa- 
tion R and will be represented by truncated power series 
in £ as the potential curves: 



a(R) 






2R„ 



R + Rr 



• E 

i=0,l,. 



Oii-C 



(8) 



U ad {R) = j 1 - — 



2R 7 , 



R + R n 



E 

i=0,l,... 



vi-e (9) 



with /z re f being the reduced mass of the selected ref- 
erence isotope, here 39 K 2 . These equations consider the 
asymptotic behaviour of the corrections, which should 
vanish for R — > oo because the intermolecular motion 
becomes infinitely slow, and for the adiabatic correction 
the proper exponent, namely 6, to keep the lowest order 
of the long range behavior as in Eq. ([5]) . The adiabatic 
corrections will be added to the potentials in Eq. (UJ 
and the operator T n will be conventionally split into the 
one dimensional kinetic energy of vibration and the ro- 
tational energy where a(R) can be added according to 
Eq. ([7]). With this we have now the complete physical 
model for calculating the necessary eigenenergies for the 
spectral data and the scattering properties to derive the 
resonances. 

The evaluation for the free parameters in the model 
is performed in a self-consistent iteration loop, in which 
the fit of the Born-Oppenheimer potentials, that describe 
the spectroscopic data with the help of the Numerov pro- 
cedure for solving the Schrodinger equation as given for 
the states X 1 S+ and a 3 £+ separately by Eq.flTJ), is 
alternated with the coupled channels calculations for de- 
scribing the Feshbach resonances and photoassociation 



data. The coupled channels calculations are also used 
for deriving the hyperfine corrections to the observed 
spectral lines and thus constructing the hyperfine free 
spectroscopic data, which are the inputs to the potential 
fits. The potentials are obtained simultaneously for both 
states, because we have data which represent the differ- 
ence between the singlet and triplet manifold of rovibra- 
tional levels and both states have a common asymptote 
with a common long range function as given by Eq. ([5]). 
This is a very important constraint for both potential 
functions, often neglected within the evaluation of sets 
of potential functions. 



V. EVALUATION AND RESULTS 

The most recent analysis of the two ground states of 
potassium is from our own group 23] . This work used the 
spectroscopic data from Fourier transform spectroscopy 
on the state X 1 T<£ from [221 ] and own measurements 
for this state and the state a 3 S+ , and furthermore, re- 
sults from Feshbach spectroscopy by [lg, LUJ 120, [2l| and 
photoassociation spectroscopy by [ll[ . With the present 
beam experiment, the remaining energy gap between 
the range covered by conventional spectroscopy results 
and those from the Feshbach resonance studies of about 
0.9 cm -1 and 1.7 cm -1 for the states X 1 E+ and a 3 E+, 
respectively, could be closed for the state X 1 S+, but no 
new measurements were obtained for the state a 3 £+ . 
However, the observed hyperfine structure of the highest 
vibrational levels in X 1 S+, see Fig. |31 is a clear signa- 
ture of the coupling between the singlet and the triplet 
state, and introduces new information. 

To learn how much could be gained by these additional 
pieces of information, we calculated the term energies of 
the newly observed level with the formerly published po- 
tential curves J23|. All new measurements showed sig- 
nificant deviations to the simulation results. First, the 
measurements contain very precisely determined energy 
intervals between v = 41 and levels from v" = 80 on, here 
the deviations were about twenty times the experimen- 
tal uncertainty and still about two times the expected 
accuracy predicted from the uncertainty of the Fourier 
transform data. Second, the differences between pairs of 
asymptotic levels, say for example v" — 80 to 83 etc., 
show deviations in the order of twenty to forty times the 
experimental uncertainty. Thus the new data give impor- 
tant information for connecting the set of deeply bound 
levels with the asymptotic ones derived from Feshbach 
resonances. 

We started the full routine of self consistent field ap- 
proach by iterating between potential and asymptotic 
fit of spectroscopic data and of Feshbach resonances, re- 
spectively. For these calculations we reduced the weight 
for the large vibrational intervals (v = 41 to v" = 80 
and above) by selecting as experimental uncertainty only 
0.001 cm -1 instead of the actual value as derivable from 
tables HUT] to be between 0.0001 and 0.0002 cm" 1 . Oth- 



erwise, those high weights might locally distort the po- 
tential fit. We already mention here, that the final fit 
result reproduces these data within the original experi- 
mental uncertainty with the exception of two cases where 
the deviations are two times the uncertainty. 

Because we are interested to use the information on 
the Feshbach resonances up to their experimental uncer- 
tainty, we checked first, if thermal averaging is needed 
to describe the observed resonance position in the ultra- 
cold ensemble at the specific conditions from the differ- 
ent laboratories. We followed the procedure described by 
Ticknor et al. [l8[ and integrated the energy dependent 
cross sections or rates of the assumed elastic collision over 
a Boltzmann distribution, because the experimental con- 
ditions were sufficiently far from the regime of quantum 
degeneracy Most profiles plotted linearly as function of 
magnetic field show an asymmetric form, but the peak 
positions deviated very little from those expected when 
using as single kinetic energy the value derived from the 
reported temperature. Thus thermal averaging during 
the fit procedure was unnecessary. We would like to men- 
tion, that the functional form of the cross section with 
energy given in Fig. lb in [18[ for s-wave scattering and 
rationalized in that paper is not always so simple. We 
found cases where also peak structures appear, namely 
if the resonance peak is located at an energy before the 
unitarity limit is approached. These cases lead to more 
asymmetry in the resonance profiles as functions of mag- 
netic field as in the cases of p-wave resonances, but we 
did not find for s-waves such severe examples as given in 
Fig. 2a of 18] for p-waves. We should remind the reader 
that we simply assume that the resonance position can 
be derived from an elastic two-body collision, whereas 
in the experiment often trap loss is observed, and thus 
a three body process is involved. Deviations have been 
studied for the case of Rb by Smirne et al. [4JJ. Ex- 
tending our analysis in this direction is far beyond our 
present capability, developments of new theoretical mod- 
eling approaches, however, would be desirable. 

We followed three different fit cases to get a reliable 
theoretical description of all observations and by com- 
paring the results we might verify the validity of the 
BOA and get a limit of its precision or derive a quan- 
titative correction to the BOA. First, we applied all data 
to fit potentials without any Born-Oppenheimer correc- 
tion. Second, we split the data set: The spectroscopic 
data and the Feshbach resonances for 39 K2 only will de- 
termine the potentials for states X 4 E+ and a 3 S+ which 
are then used to simulate the few Feshbach resonances 
for 40 K2. Our attention was concentrated on any sys- 
tematic shift seen between the observed and simulated 
resonances for 40 K2. Third, we included correction func- 
tions as shown in Eq. §§§. We again include the 40 K2 
data and check if any improvement on the fit compared 
to the first case can be obtained. 

For the first case the fit of the whole data set starts 
with the potential approach reported in [23j |. the num- 
ber of potential parameters for the intermediate regime, 



TABLE IV: Parameters of the analytic representation of the 
X 1 Sj state potential without any Born-Oppenheimer cor- 
rection. The energy reference is the dissociation asymptote. 
Parameters with * are set for continuous extrapolation of the 
potential. 
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Eq. ([3]), and the values of -Rinn and i? ut for splitting the 
R axis into the different regions are taken from that ref- 
erence. The new measurements for v" > 80 have outer 
turning points in the region beyond i? ou t thus give ad- 
ditional information on the long range form. This is 
especially important because several different rotational 



TABLE V: Parameters of the analytic representation of the 
a 3 Y, U ~ state potential without any Born-Oppenheimer cor- 
rection. The energy reference is the dissociation asymptote. 
Parameters with * are set for continuous extrapolation of the 
potential. 
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states are observed. This gives a more detailed test of the 
potential function than the vibrational spacing of levels 
with the same rotational angular momentum alone. For 
the potential fit, the calculated hyperfine structure of 
these observed levels is subtracted, which was derived in 
the preceding iteration step with coupled channels calcu- 
lation. Typically, three iteration loops, alternating po- 
tential fit and coupled channels calculation, are needed 
to get convergence. 

Tables IIVI and [V] show the fit results of the first case 
without any attempt of Born-Oppcnhcimcr correction for 
the two states X *E+ and a 3 E+. The total fit of the 
Born-Oppcnhcimcr potentials contains 8775 data points, 
in which Feshbach resonances are included by convert- 
ing them to extrapolated rovibrational levels of those 
quantum states which produce the observed resonances. 
These are vibrational levels with v£ = 84 and 85 for 39 K2 



and 40 K2, respectively and u" = 26 for both isotopomers. 
Because for 40 K2 a p-wave resonance was also reported, 
a J" = 1 and A^" = 1 level for those vibrational states 
of X 4 E+ and a 3 E+ , respectively, was constructed. A 
normalized standard deviation of a — 0.78 was obtained 
for the potential fit in the iteration loop. 

In the last iteration step the fit of the 12 Feshbach res- 
onances yields a standard deviation of a = 0.76, which is 
very satisfactory from the statistical point of view, but 
the detailed inspection of the distribution of deviations 
between the two isotopomers shows that for the four res- 
onances of 40 K2 a smaller scatter of the deviations is ob- 
tained than in the case of the resonances of 39 K2. This 
asymmetry indicates that the quality of the fit might be 
not as good as the single number of a would tell. Hence, 
we have good reason to try the second case of our fit 
strategy. 

But before describing these results we would like to 
mention, that the photoassociation measurements on 
39 K2 from Ref. [TlJ were incorporated in the fit and 
showed good consistency with the other data. On av- 
erage, these measurements are at least a factor of 10 less 
accurate than the Feshbach data, if one transforms by the 
Zeeman effect the magnetic field accuracy to an uncer- 
tainty of level energies. Thus these measurements have 
a small weight in the evaluation. For all later calcula- 
tions below, the photoassociation results are always in 
good agreement, thus we do not need to mention it later 
again. 

In the second case we fitted only the isotopomer 39 K2 
and used the result for calculating the four Feshbach res- 
onances of 40 K2. This case corresponds to the scenario of 
many cold collision studies, where Feshbach resonances 
are observed for one isotope combination and mass scal- 
ing is applied for predicting resonances of other isotope 
combinations. From the statistical point of view the over- 
all fit quality in this second case was as good as in the 
former one. The resonance fit of 39 K2 is of particular 
interest which yields a a — 0.55. This small value ob- 
tained by restricting the data set to a single isotopomer 
directly reflects the above mentioned asymmetry in the 
distribution of the two isotopomers. With these poten- 
tial functions for X 1 Ej" and a 3 E+ the simulation of the 
resonances for 40 K2 results in deviations of more than 
0.1 mT and all are negative. Because the experimen- 
tal uncertainties of these resonances are mostly below 
0.01 mT this second fit cannot be accepted for the de- 
scription of isotopomer K 2 from results of isotopomer 
39 K 2 . 

This gives us good reason to turn to our third case, in 
which we introduce Born-Oppenheimer correction func- 
tions. The very sensitive point when introducing new 
functions is the decision how many terms, represented 
by the number of additional parameters, should be tried. 
Indeed, we are looking for small effects buried in a large 
set of data, which is already described fairly well in the 
first case, and we have to expect correlation between the 
potential functions and new correction functions. Thus 



10 



we concentrated only on the adiabatic correction, [/ a( j 
of Eq. $)§, with two expansion coefficients for the state 
X ^ij" and a single one for the other, because here the 
isotope variation is only obtained from the resonances 
and thus for the asymptotic levels. We selected 39 K2 as 
the reference isotopomer, because this species had the 
largest data set. 

The potential fit yields the same quality as for the 
former ones, namely a standard deviation a = 0.78, but 
the fit of the Feshbach resonance of both isotopomers 
results in an improvement with a standard deviation of 
a = 0.59 compared to the value 0.76 obtained in the first 
fit. No asymmetry is found in the present distribution 
of the deviations. Tables IVII and IVIII show the new fit 
results. The correction function for X 1 T,+ is a linear 
function in the expansion variable £, compare Eq. ^. 
We also tried a fit by using a quadratic function with 
two parameters, i.e. with Do x £° and V2 X £ 2 , which will 
change the influence mainly for the low levels within the 
singlet potential. But here the potential fit was not as 
good as before, which indicates that the few spectral lines 
to deeply bound levels of the singlet ground state for the 
isotopic species 39 K 41 K also play a non-negligible role in 
the final evaluation. 

The influence of the derived quantitative Born- 
Oppenhcimer correction can be seen by comparing the 
term energies of the obtained model to the term ener- 
gies obtained from the same model with the corrections 
artificially set to zero. Figure 2] shows these energy dif- 
ferences of the vibrational ladders for 40 K2. A graph 
for 39 K 41 K would almost coincide with that shown be- 
cause the reduced masses of these two isotopomers are 
very close in value. The differences of the vibrational 
ladders are proportional to the expectation values of U a d 
for the vibrational states v of the state X 1 S+. For low 
vibrational levels the differences are on the order of the 
experimental accuracy. For the asymptotic vibrational 
level they almost vanish as desired by fixing the asymp- 
totic behavior according to the long range function in 
lowest order in Eq. ©. The actual functional form seen 
in Fig. 3] is far from being fixed by this study, only 
the magnitudes should be considered. For v" — 85, the 
level which correlates to observed Feshbach resonances in 
40 K2, the energy difference is about 2 MHz and thus at 
least a factor 3 larger than the uncertainty of the reso- 
nance determination would tolerate and it corresponds to 
the resonance shifts for 40 K2 between the observed val- 
ues and those calculated with the potentials derived form 
for 39 K2 data as obtained in the second fitting case. A 
similar magnitude results for the levels of state a 3 £„ . In 
total, we can state that we reached the limit of accuracy 
of the Born-Oppenheimer approximation in modeling the 
observations on K2. Small improvements of the experi- 
ments will allow to derive clearly the magnitude of the 
corrections to the Born-Oppenheimer approximation. 

Finally, we also checked if retardation effects could al- 
ter such statement about accuracy and possible correc- 
tion to the Born-Oppenheimer approximation. As de- 
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FIG. 4: Difference of vibrational ladders for K2 calculated 
with the potential from table IVTl with adiabatic correction or 
setting it to zero. 



scribed by Marinescu et al. [42J we introduced the cor- 
rection factors calculated by them for all C{ coefficients 
in Eq. $5§ and repeated the iterative fit. We got the 
same quality of the fit and all changes were insignifi- 
cant regarding the conclusion on the magnitude of possi- 
ble Born-Oppenheimer corrections and derived scattering 
properties, which will be collected in the next section. 



VI. 



CONCLUSION 



The new measurements with a molecular beam setup 
close the level gap between the data from Fourier trans- 
form spectroscopy and Feshbach resonances. The simu- 
lations of the new levels with the earlier derived poten- 
tials [23| showed that this interpolation results in pre- 
dictions, which deviate by about two times the error 
limit from the Fourier transform spectroscopy. Precise 
measurements by Feshbach spectroscopy do not help sig- 
nificantly for obtaining better predictions in the range 
below those vibrational levels, which are determined by 
the resonances. This is important information for other 
studies in systems, where a combined analysis of data 
sets from different, not linked energy regions is per- 
formed. A direct comparison between the earlier poten- 
tial curves and the present ones gives differences on the 
order of few times the uncertainty of the Fourier trans- 
form data, more specificly about 0.02 cm -1 and for the 
steep repulsive branch at short internulcear separation 
more, distributed over the full potential curve. Solv- 
ing the Schrodinger equation means integration over the 
full range R of the potential function and results for the 
eigenvalues only in changes of about 0.003 cm -1 or less 
as seen between the earlier simulated asymptotic level 
energies and those measured now. One should also keep 
in mind that the A-scheme of excitation applied for these 
new measurements gives much improved accuracy at an 
isolated deeply bound energy regime, i.e. around v = 41. 
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This fact certainly influences also the overall change in 
the potential function. We conclude, that having data 
for the full range of potential energy with no gaps will 
give a secure basis for reliable predictions, whereas ex- 
trapolations and even interpolations should be labeled 
by a remark of caution for using those. Our analysis also 
shows how large changes in potential energies can ap- 
pear compared to the overall accuracy of level energies, 
namely in our case about a factor of ten. 

Such limitations also restrict the chances to discuss 
higher order corrections in the theoretical models like 
our desire for searching for Born-Oppenheimer correc- 
tions. We think that we are here in a favorable situa- 
tion because at least for one isotopomer, namely 39 K2, 
we have spectral information over the full range of the 
potential energy, including the atomic asymptote by the 
Feshbach spectroscopy. Comparing the different fit re- 
sults from the three cases, (i) evaluation with no cor- 
rection, (ii) separate evaluation of the isotopomers and 
finally (iii) a simultaneous evaluation including Born- 
Oppenheimer corrections, we found the first indications 
of Born-Oppenheimer corrections for the bound levels 
correlating to the Feshbach resonances. 

For characterizing the cold collision properties at zero 
kinetic energy we calculated the scattering lengths of 
the singlet and the triplet state from the results of the 
first and third case of fitting for the isotopic species of 
diatomic potassium of natural abundance. Table IVIIII 
shows the results and includes also reports from other 
authors. The error limits are inferred by conventional 
methods including the covariance matrix from the non- 
linear least squares fit in the linear limit; they should be 
interpreted as ler standard deviation from the statisti- 
cal distribution of the deviations within the fit. For the 
comparison of the different approaches we discuss only 
the results for 39 K 2 and 40 K 2 , because in these cases we 
have direct experimental observations. All others are de- 
rived quantities only, and the respective model was used 
for their calculations. For the singlet state the scattering 
length of 39 K2 is different slightly outside the statistical 
uncertainty limit between the approach with and without 
Born-Oppenheimer correction. For the triplet state the 
same appears for the isotopomer 40 K2. Because one sees 
this difference only slightly outside the statistical uncer- 
tainty limit we cannot take this information as a conclu- 
sive observation of a correction to the Born-Oppenheimer 
approximation for the modeling of cold collisions, but we 
believe that it gives the limit of mass scaling and that 
any improvement in precision would need to include cor- 
rections in the theoretical description like in the form of 
Eq. ([7]) . Cold collisions outside resonance cases are prob- 
ably less influenced by possible corrections to the Born- 
Oppenheimer approximation and therefore mass scaling 



of scattering lengths might be widely applicable, but res- 
onances themselves seem to be more sensitive, because 
we observed that mass scaling from 39 K2 to 40 Fv2 shifts 
resonances in the order of 0.1 mT, hence outside the most 
recent experimental error limits. More high precision ex- 
periments are desirable to obtain a definite conclusion. 

The results on scattering len gths are mostly in good 
agreement with our earlier ones J23il , small deviations to 
those reported by D'Errico et al [2l| could stem from the 
different fit qualities as already noted in [23j |. 

New precise studies of different isotopomers in other 
molecules are underway. We learned from Simoni [44] 
about new measurements on 39 K Rb, which they com- 
bined with earlier results on 40 K 87 Rb [H H| HI . They 
concluded, that the resonances can be described apply- 
ing mass- independent potentials for both states X 1 E + 
and a 3 E + . But the derived potentials concentrate only 
on the asymptotic behavior to get a predicting model 
for ultracold collisions and the capability to precisely re- 
produce more deeply bound levels is lost. In this way 
one obtains two solutions with different ranges of valid- 
ity, namely the deeply bound range [35| or the asymp- 
totic range. We believe that it would be valuable to re- 
peat the procedure like in the present study of K2 with 
full potentials and to keep the whole data set from spec- 
troscopy and Feshbach measurements simultaneously in 
the fit. Otherwise, small corrections can be overseen in 
the artificial separation of the two solutions. 

Similarly, new studies for LiK evolve. Feshbach spec- 
troscopy was reported for the collision 6 Li and K by 
the Innsbruck group [6| and Fourier transform spec- 
troscopy was continued in a Ph.D. thesis by Houssam 
Salami [48( giving a large amount of data for 6 Li 39 K 
and 7 Li 39 K for which Salami derived corrections of the 
Born-Oppenheimer approximation for the deeply bound 
region of X 1 S + . Here we have the favorable case that 
the reduced mass of the system is mainly determined by 
the light atom Li and thus the mass effect from 6 Li to 
7 Li is larger than for heavier systems. First steps done 
by us with our approach already indicate very promis- 
ing results, which could show the variation of the Born- 
Oppenheimer correction from small internuclear separa- 
tions to the asymptotic range, something that could only 
be guessed from Fig. 0] in the case of K2. 
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TABLE VI: Parameters of the analytic representation of the 
X E+ state potential with adiabatic Born-Oppenheimer cor- 
rection and reference isotopomer 39 K2. The energy reference 
is the dissociation asymptote. Parameters with * are set for 
continuous extrapolation of the potential. 



R < Rinn = 2.870 A 


A* 


-0.263145571 


x 10 4 


cm 1 


B* 
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-1 A 12 
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A 
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Rm 


3.92436437 A 
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0() 


-4450.899484 


cm -1 


«1 


0.30601009538111 x 
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a 3 
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TABLE VII: Parameters of the analytic representation of the 
a 3 S+ state potential with adiabatic Born-Oppenheimer cor- 
rection and reference isotopomer 39 K2. The energy reference 
is the dissociation asymptote. Parameters with * are set for 
continuous extrapolation of the potential. 
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TABLE VIII: Scattering lengths (unit a = 0.5292 A) for the 
both approaches of the potentials (A: no B.-O. correction, B: 
with B.-O. correction) for different isotopomers of potassium. 
Results from experimental analysis by other authors are given: 
a EL bl. 
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